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Abstract 

We study the nucleon at finite temperature using Bethe-Salpeter equation (BSE) where it is 
considered as a bound state of a scalar-diquark and a quark. The interaction between the diquark 
and the quark is taken as an exchange quark. The constituent quark and diquark substructures are 
considered momentum and temperature dependent based on the lattice QCD and dynamical chiral 
symmetry breaking (DCSB). The kernel interaction modification with temperature is approximated 
using the mixed representation in the imaginary time formalism and then we adopt an adiabatic 
(instantaneous) approximation after evaluating the Matsubara frequencies sum. The results of 
BSE are compared with those for MIT bag model. We argue the bag pressure is a temperature 
dependent and it decreases smoothly with temperature to be consistent with the results for BSE. 
Furthermore, the nucleon bound state obtained by BSE ceases to exist just below the critical 
temperature and this indicates that the nucleon deconfines to its constituent quarks at the critical 
temperature. The nucleon mass decreases with temperature as Mn(T)/Mn(0) ~ 1Z S (T) where 
TZ(T) = (l - b^ 1/2 and 5 « 1 and this suggests B(T)/B(0) = TZ 45 {T) for MIT bag. 



I. INTRODUCTION 



The properties of hadrons and their masses at finite density and temperature are ex- 
tremely important to study the phase transition from the hadronic phase to the quark gluon 
plasma. Generally speaking, the nonperturbative calculations at finite temperature and/or 
density support the idea that hadronic masses change as a function of temperature and 
density. The nucleon mass is found to depend substantially on temperature variations of 
the quark condensate and interaction process. As the condensate disappears smoothly with 
chiral-symmetry restoration so does the nucleon mass 

Recently, baryons have been studied in the context of Bethe-Salpeter equation (BSE). 
The baryon is considered as a bound state of the confined constituent diquark and quark 
interacting via an exchange quark. Despite the simplicity of the model, it has been used 
to study the intermediate energy physics successfully Furthermore, it has 

been some attempts to extend BSE to study mesons at finite temperature jlO . 

HQ- How- 
ever, the dressed quark propagator and the variation of the constituent quark mass with 

temperature are essential to study the bound states in a hot bath. The Schwinger-Dyson 
equation (DSE) is used to study the quark propagator at finite temperature The DSE 
with QCD potential can be solved numerically at finite temperature Q|. Furthermore, the 
constituent quark mass and its variation with temperature is studied using the lattice-QCD- 
based Schwinger-Dyson equation. Spontaneously broken chiral symmetry is found to be 
restored at the critical temperature [3, Q]- The linear confinement dissociates smoothly 
with temperature From DSE, it is expected that the dressed quark mass to decrease 
with temperature in the same order of the linear confinement dissociation. The condensate 
is found also to disappear smoothly at the same order. 

In Sec. II, we shall derive the BSE for a quark-diquark bound state at finite temperature. 
The modification of the interaction with temperature is approximated using a mixed repre- 
sentation of the imaginary time formalism and then we adopt the adiabatic (instantaneous) 
approximation. In this approximation, we set the time component to zero after evaluating 
the Matsubara sums not before. The constituent scalar-diquark and quark substructures 
are effectively taken into account and their masses are taken momentum and temperature 
dependent based on the lattice-QCD. In Sec. Ill, we briefly review the MIT bag model for 
nucleon and then we extend it to finite temperature case by taking on the considerations 
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the dissociation of the nucleoli bag pressure with temperature. Finally in Sec. V, we present 
our results and conclusions. 



II. DIQUARK-QUARK BETHE-SALPETER FOR NUCLEON 

The baryon is considered as a bound state of confined constituent diquark and quark 
interacting via quark exchange. The BSE for the diquark and quark bound states is given 
by 

d 4 p' 



where 



Xp(p) = ^D(Pi)S q (p 2 ) 



(2tt)< 



[-iK(P,p,p')} X p( P '), 



(1) 



Pl = TjiP + p 
p 2 = 1] 2 P ~ p. 



(2) 



')S q (q) where <fi(q) is the quark-diquark vertex coupling 



We have taken K(P,p,p') 
gj) q = g^ q (q) and S q (q) is the exchange quark propagator and q = p + p'. The scalar diquark 
propagator reads 

A D (pi) = 



Z D J \pj-mj + i5} 
Z D J 2ei 



Po + M 1 -e 1 +i5 po + Mi + Cl - iS 
and the fermion quark propagator reads 

1 \ 7-^2+^2 



(3) 



SM = 



—I 



Z q ) p\ — ml + iS 

M-p) 



+ 



A-(-p) 



Po - M 2 + e 2 - iS po - M 2 - e 2 + id 



(4) 



where q = y p 2 + mf and i — 1 refers to the scalar diquark while i — 2 refers to the quark 
and Mi = r^M and M 2 = (1 —r)i)M. We have chosen the proper choice rji = mi/{mi+m 2 ). 
In the center of mass frame we have P = (M, 0). The fermion positive and negative energy 
projections read 

e 2 + a/3(— 7 • p + m 2 ) 



M-p) 



2e 2 



(5) 
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where a = ±. It is possible to take Zd = 1 and Z q = 1 with a good approximation. In 
principle, these functions can be considered in the calculations by adopting smearing func- 
tions in the propagators j?], Q]. In the present calculations, we shall assume these functions 
are absorbed by the quark-diquark vertex coupling g 2 Dq . The constituent diquark and quark 
masses are given by mjj = m^p) and m q = m g (p) for the diquark and quark, respectively. 
In the adiabatic approximation, the temperature response is separated from the momentum 
dependent and the constituent diquark and quark masses reduce to = m£>(p;T) and 
m q = m q (p;T). The Dynamical chiral-symmetry breaking (DCSB) is one of the most non- 
perturbative QCD feature. Based on the Lattice QCD at the quenched level and Landau 
gauge and the lattice data in the chiral limit, the quark mass function at zero temperature 
is reproduced by 

muAp ) = mo+ [i+(pi/p 2 )^t (6) 

where p = 870 MeV, 7 = 3.04, the current mass m = MeV and p\ — p\ + p 2 is the 
momentum in Euclidean space. Using the lattice- Q CD-based gap equation (DSE) at finite 
temperature, it is found that rriQ(T) has a chiral symmetry restoration phase transition at 
the critical temperature T c . The quark and diquark acquire their constituent masses from the 
linear confinement. These masses can be computed using DSE to calculate the self-energy 
correction. However the linear confinement constant V(r) ~ a^^r in the configuration 
space is found to dissociate thermally spontaneously to ar (T)/clt (0) = H(T). The linear 
confining is a nonperturbative effect and it is believed to generate the major contribution 
to the self-energy correction for the current quark mass to acquire its constituent value. 
It is reasonable to take the self-energy correction and mass variation with respect to the 
temperature at the same order of the spontaneous linear confinement dissociation 

m Q (T)/m Q (0)=TZ(T), (7) 

where K{T) = y/l - bT 2 /T? for T < T c and K{T) = \fl^b for T > T C Q, 20]. Since the 
constituent diquark is assumed as a weekly quark-quark bound state, it is reasonable to as- 
sume mj3 (T)/mj}(0) = TZ(T). Furthermore, we assume the diquark substructure takes the 
same quark formula substructure given in Eq.fjHJ). In the adiabatic approximation, euclidean 
momentum p\ — p\ + p 2 in the quark substructure is reduced to three dimensional momen- 
tum where p\ = 0. However, in the imaginary time formalism, we can write p^ = (2n + l)7rT 



where n — 1, 2, • • • and we can take p\ = ir 2 T 2 + p 2 to the first order approximation. To be 
consistent with the adiabatic approximation considered in our calculations, we have adopted 
pi = in our calculations. Furthermore, we have shown that the first order approximation 
effect is small. 

In an adiabatic approximation, the bound state equation at finite temperature is derived 
as follows 

d 3 p' 



dpo 

(27T)' 



dp'p 

(2tt)' 



(8) 



At zero temperature, the adiabatic approximation is reduced to the instantaneous approxi- 
mation. However, as we shall show below, a special attention must be paid for the interaction 
potential in the adiabatic approximation. The kernel interaction is taken as an exchange 
quark between the diquark and quark, 



-iK(P;p,p')xp(p') = -iK(q)xp(p') = g 2 Dq S q (p + p')xp(p), 

where q = p + p'. The BSE propagator is reduced to 

dp 



£(M,p;T) = i J J^A D ( Pl )S q (p 2 ) = 

= g + (M,p;T)A + (-p)P + g_(M,p;T)A- 



(9) 



(10) 



where 



£ + (M,p;T) = i- 



1 



+ 



2e 1 



M - (ei + e 2 ) 
1 

M+{e 1 - e 2 ) 



[1 +n B (e_i) - n F (e_ 2 )] 
[riB (e+i) + n F (e_ 2 )] 



(11) 



and 



(?_(M,p;T) = — 
1 



M+ (e 1 + e 2 ) 
1 



[1 + n B (e+i) - n F (e +2 )] 
[«fl (e-i) +n F (e+2)] , 



(12) 



M-(ei-e 2 )_ 

where e±j = ± Mj. The thermal distribution functions for quark and scalar diquark 
n F {x) = eX /T +1 and ub(x) = eZ/T _ 1 are the Fermi-Einstein and Bose- Einstein distribution 
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functions, respectively, and x > 0. At zero temperature, the BSE propagator given by 
Eq. ffTUj) is reduced to 



^A 1 ( Pi )S 2 (p 2 )\t=o = tt 

Z7T Z6i 



+ 



(13) 



_+M-(ei + e 2 ) +M+(e 1 + e 2 ) 
It is a wonderful to write the thermal distribution function of BSE propagator in a proper 
way to discuss thermal reaction processes. The thermal distribution function for the positive 
energy solution is written as follows 



[1 + n B (e-i) - n F (e- 2 )] = [1 + n B {e-\)) [1 - n F (e- 2 )) + n B (e^ 1 )n F (e- 2 ) 



(14) 



It corresponds to the difference between the virtual decay and creation rates via the virtual 
reaction processes 4>m — ► 0102 and (fi±(fi 2 — > 0a/, respectively It is interested to note that 
the constituent diquark and quark chemical potentials Mi and M 2 , respectively, are always 
less than their relativistic energies and e_i = e% — Mi > and t- 2 = e 2 — M 2 > to preserve 
the confinement condition. This guarantees the creation rate is always larger than the decay 
rate. Furthermore, the additional 'Landau Damping' term 



[n B (e +1 ) + n F (e_ 2 )] = n B (e + i) [1 - n F (e^ 2 )} + [1 + n B {e +1 )\ n F (e- 2 ) 



(15) 



corresponds the virtual decay and creation processes via the reaction scattering 4>i4>m ~~ * 02 
and 02 - > 0i 4>m, respectively. This corresponds particles disappear or are created through 
the scattering in the bath, and not via the process which are available at zero temperature. 
The thermal distribution function for the negative energy solution corresponds to <p M 



'1^2 



and 



for the decay and creation processes, respectively. On the other hand, the 



'Landau Damping' term in the negative energy solution reduces to the creation and decay 
processes 2 — > 0i0a/ and 4>i4>m ~ > 02; respectively 

The exchange quark propagator is decomposed as follows 

A + (-q) , A_(-q)' 



S q (q) = ifi 
Hence the kernel interaction reads 



- K(q , q) = g 2 D J3 



A, 



qo 



qo + eq 



lS) + A-( 



-qj 



Qo 



Z Q J 



When we substitute q = 0, we get 
-K{q o = 0,q) 



-9 2 Dq P 



A- 



-qj 



A_( 



-qj 



(16) 



(17) 



(18) 
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However, this choice is not adequate in the adiabatic approximation. In the adiabatic as- 
sumption, we shall consider the instantaneous approximation after evaluating the Matsubara 
frequencies sum in the mixed representation for the imaginary time formalism. The Fourier 
transformation of the kernel interaction to the mixed representation reduces to 



K(r,q;T) = T ^ e~ g ° T K(qo, q) 



= ~9d. 



q [e^ T [1 - n F (e q )} /3A + (-q) + [n F (e q )} /3A_(-q) 



(19) 



(20) 



The inverse transformation of the mixed representation with an adiabatic approximation 
reads 



K(q Q , q; T)| ?0=0 = f T dre^\ qo=0 K(r, q; T) 

J 



= -g 2 D - [1 - 2n F (e q )} (5 [A+(-q) + A_(-q)] , 
e q 

= -9d q — t 1 - 2n F (e q )]P, 



(21) 



where e q = ^/(p + p') 2 + "^ 2 - It is worth to note here that in the adiabatic approximation 
we have substituted qo = in the mixed representation after evaluating the Matsubara 
frequencies sum. We have introduced in BSE the following interaction potential 



-K(q ,q;T)\ qo=0 = V(q;T)P, 
where V(q; T) = -g 2 Dq j- q [1 - 2n F (e,)]. The BSE reduces to 

A<Kp)U =1 = S(M,p;T)/3 / |^(q;T)0(p'). 

The BSE wavefunction decomposes as Dirac wavefunction for positive parity 

/ 



(22) 



(23) 



The BSE wavefunction reads 



0(P) 



( 



\ 



Vij{P) 



V^-pOj(p) 



jrn' 



(24) 



+ (p) 

r(p) 



v o- • pC(p) 



X, 



(25) 
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where x is the spinor and the quantum numbers are suppressed. The decompositions of the 
BSE wave function with respect to the positive and negative energy components read 

/ 



A + (- P )0(p') = ^ 
2e 2 



[e 2 + m 2 ]v(p') ~ |p|[p p' + ipxp' ■ a]((p') 



\<7-P 



IpI^(p') + h - m 2 ] [p ■ p' + ip x p' • 5]C(p'; 



and 



X, (26) 



A-(-p)0(p) 



1 

2T 2 



[e 2 - m 2 ]r/(p') + |p| [p-p' + ipxp'- <?]C(p') 



\^-P 



IpI^(p') + h + ma] [p • p' + ip x p' • a]((p') 



(27) 



We have taken g 2 D = g 2 , e A2q2 where q = p + p' and A = 0.18fm 2 , g dqq = 14.150], 



mi — 650 MeV and m 2 = 350 MeV for the diquark and quark, respectively. 



III. MIT BAG MODEL FOR NUCLEON 



The nucleon is assumed as MIT bag of the confined quarks 
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□ 



17, 18]. The quark field 



if) q (r,t) inside the nucleon bag of radius R satisfies the Dirac equation 



i^d^ -mjl ip q (r,t) = 0, 



(28) 



where m° q is the current mass of a quark of flavor q. The current quark masses are taken 
m u = rrid = for the up and down flavor quarks. The single-particle quark energy is given 
by 



(29) 



where 



x q + R?m 2 . 



(30) 



For a given value of the bag radius R, the quark momentum x q is determined by the confine- 
ment boundary condition at the bag surface which, for quarks of flavor q in a spherical bag, 
reduces to jo(x q ) = (3 q ji(x q ), where (3 q - 



n q +Rm an d io> ji are spherical Bessel functions 
of order zero and one, respectively. The bag energy is given by 



Z 47T 3 



(31) 
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where Ylq q -r is the total constituent quarks kinetic energy inside the bag and ^ term is 
the zero-point energy for quarks and B is the bag parameter. The spurious center-of-mass 
energy is subtracted to obtain nucleon mass, 



M * = V^bag - < Plm >, (32) 

where 

n q 

<Plm>=J24/ R2 - (33) 
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The bag radius R is obtained through the minimization of the baryon mass with respect to 
its bag radius 

The bag parameter B is taken as B corresponding to its value for a free nucleon. The values 
B^ 4 = 188.1 MeV and Zq = 2.03 are chosen to reproduce the free nucleon mass Mn at its 
experimental value 939 MeV and bag radius Rq = 0.60 fm. The current quark mass m q is 
taken equal to zero. 

In order to consider the thermal effects for the nucleon mass M(T), we assume that the 
bag parameter dissolves with temperature as follows 

B = B(T) = B x TZ 4S {T) (35) 

where 

K(T) = (l - b— J . (36) 

With a choice 5=1, the nucleon mass at finite temperature reduces to 

M N (T) = M N (0) x K{T). (37) 

This choice is relevant to the thermal dissociation of the linear confinement cr(T) = (TqH(T) 
[19! and its possible connection to the bag pressure dissociation B^T) = 5 1/4 ^(T). 
This connection leads to the choice 5=1. 
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IV. RESULTS AND CONCLUSIONS 



We have studied the nucleoli as a bound state of a scalar-diquark and a quark by diag- 
onalizing the BSE at finite temperature. The interaction is taken as a quark exchange in 
the adiabatic approximation. The quark and scalar-diquark substructures and their modifi- 
cations with momentum and temperature are introduced self-consistently in the BSE. The 
interaction modification with temperature is approximated using a mixed representation of 
the imaginary time formalism in the propagator and then the adiabatic (instantaneous) ap- 
proximation is adopted by substituting the time component to zero q = after evaluating 
the Matsubara frequency sums not before. We have included explicitly in our calculations 
the smooth decreasing of the constituent quark and diquark masses with temperature. The 
thermal variations of the constituent particle masses are based on the lattice-QCD calcula- 
tion. In Fig. (I), we display the dependence of constituent quark and scalar-diquark masses 
on their substructure momenta at zero temperature. The quark and diquark acquire their 
constituent masses 350 and 650 MeV, respectively, at zero momentum. They decrease with 
momentum to reach their current masses at large momentum. The diquark acts as a weakly 
bound state of quark-quark and its mass varies with temperature in the same order of varia- 
tion of its constituent quarks. In Fig. (II), we display the quark and diquark masses versus 
temperature at zero momentum. It is shown that the constituent masses decrease with 
temperature as 1Z(T). The thick lines correspond p± = ttT, while the thin lines correspond 
P4 = 0. It is shown that the first order Matsubara frequency correction has a very small 
correction to the adiabatic approximation. 

We have calculated the nucleon mass by diagonalizing BSE with several temperatures. 
The nucleon mass versus temperature is displayed in Fig. (III). We have found that 2M N (T) 
is always less than 3m^ g (T) until the temperature reaches just below the critical one. This 
means that the nucleon bound state persists to exist until the temperature reaches a critical 
one. Therefore, the quarks continue to confine in the hadronic phase and they shall not 
deconfine spontaneously until the temperature reaches the critical one. Furthermore, we 
have found just below the critical temperature, the nucleon bound state ceases to exist. It 
deconfines to its constituents spontaneously at the critical temperature. In our calculations, 
the critical temperature is taken T c = 170 MeV. This probably a signature for the first order 
phase transition. Furthermore, the nucleon mass decreases smoothly with temperature in 
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FIG. 1: The constituent quark and scalar-diquark masses versus momentum at zero temperature. 

the same manner the constituent quark mass decreases with temperature and at the same 
order of the linear confinement dissociation in the lattice-QCD calculation. These results are 
compared with those for MIT bag model. We assume that MIT bag pressure is temperature 
dependent and it dissociates as S 1 / 4 (T) = B^ 4 1Z(T) and 5 = 1. This dissociation is at the 
same order of the linear confinement dissociation in the potential models and lattice-QCD 



calculations 
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20]. We have found with this choice, the nucleon mass decreases smoothly 
with temperature as Mn(T)/Mn(0) = TZ(T). Furthermore, the choice 5 = 1/4 doesn't 
match the mass spectrum for BSE. 

Furthermore, the solutions of BSE and MIT bag are in a good agreement with each other 
but the MIT bag model lacks the first order phase transition that found for BSE in the 
critical temperature. However, the MIT bags might have the first order phase transition 
when they are embedded in a nuclear matter because of their interaction with the medium 
and the the production of particles and antiparticles. Finally, in our conclusions, we argue 
the nucleon mass in a hot bath should be a temperature dependent and it decreases smoothly 



as M N (T) = M N TZ(T) where 1Z(T) = (l - bff' 2 . 
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FIG. 2: The constituent quark and scalar-diquark masses with p = versus temperature. 




FIG. 3: Nucleon mass calculated using BSE and the scalar-diquark constituent mass at p = 
versus temperature. It is shown that the bound state persists to exist until the temperature reaches 
just below the critical one and then the bound state ceases to exist at the critical one. 
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FIG. 4: The nucleoli mass Mjv(T) versus temperature. The nucleon mass is obtained by BSE and 
MIT bag models with several sets as explained in the text. 
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